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Abstract 


Aerodynamic coefficients are presented for perfect-gas and equilibrium-air 
solutions of the Navier-Stokes equations about the Stardust reentry-capsule fore- 
body at Mach numbers of 4.6, 7, 8.5, and 10. A comparison with Newtonian- 
flow assumptions indicates a divergence of the aerodynamic coefficients from 
Newtonian-flow for Mach numbers less than 10. The static stability of the fore- 
body is reduced by a factor of 2.5 with decreasing freestream Mach number 
between Mach 10 and 7. 


Nomenclature 

A Area, m 2 

C Force or moment coefficient 

Cp Pressure coefficient 

c Reference length, m 

e Unit vector, m 

h Altitude, km 

i,j Computational indicies 

M Mach number 

n Unit normal vector, m 

P Pressure, Pa 

5 Reference area, m 2 

T Temperature, K 

t Time during trajectory, s 

t Tangent vector, m 

V Velocity, m/s 

x,y,z Cartesian body axes, m 

a Angle of attack, degrees 

7 Ratio of specific heats 

<p,0 Surface inclination angles, rad 
p Density, kg/m 3 

Subscripts: 

.4 Axial force 

E Surface element 

N Normal force 

M Pitching moment 

O 2 Post-shock stagnation value 

00 Freestream 


Introduction 

Stardust [1] is a comet sample-and-return mission slated for the end of the 
century. Trajectory calculations for Earth atmospheric reentry require a deter- 
mination of the capsule aerodynamics. Of primary interest for the axisymmetric 
capsule are lift, drag, and pitching moment coefficients. 
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Laminar, thin-layer Navier-Stokes solutions are obtained at four trajectory 
points, Mach 4.6, 7, 8.5, and 10, for the Stardust forebody, assuming a base pres- 
sure equal to freestream static pressure. Use of this base pressure assumption 
at similar Mach-number conditions for the Commercial Experiment Transporter 
(COMET) reentry capsule has been shown to produce errors of one percent or 
less[2, 3]. A complete set of Stardust aerodynamics has been published by 
Mitcheltree et a/[ 4], including a subset of the present data. 

Navier-Stokes solutions are compared with results using Newtonian-flow as- 
sumptions, revealing a sudden shift from Newtonian-like flow at the highest 
altitude/Mach number point to a non-Newtonian surface pressure distribution 
at lower altitudes/Mach numbers. While the trend away from Newtonian-flow 
as the Mach number decreases is expected, of note for this case is how abrupt 
the departure is between Mach 7 and 10. Over this Mach number range, cov- 
ering eight seconds of flight time, the forebody static stability decreases by a 
factor of 2.5. 


Configuration 

The physical configuration of the Stardust reentry-capsule forebody is a 
spherically-blunted 60 degree half-angle cone. The nose radius is 0.23 m, the 
shoulder radius is 0.02 in, and the base radius is 0.41 m. 

The axisymmetric computational grid contains 30 streamwise cells and 64 
cells normal to the body, in a structured quadrilateral framework. Grid adaption 
was performed to align with the bow shock and boundary layer. 

For the three-dimensional cases, a singularity-free grid is employed with a 
58 by 29 cell surface mesh and 48 cells normal to the body. Figure 1 displays 
a side view of the surface mesh. Only the starboard side of the vehicle was 
computed, with left-right symmetry assumed in the flowfield. Adaption v/as 
performed on this grid to align with the bow shock and boundary layer. 

Computational Methods 

Laminar, thin-layer Navier-Stokes calculations were performed using the 
Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA) [5, 6]. 
LALTRA is a second-order accurate finite- volume approximate Riemann solver [7]. 
Previous applications to Earth reentry capsules include COMET[2, 3] and the 
Aeroassist Flight Experiment [8]. Both perfect-gas and equilibrium-air[9] models 
were employed. 

Modified-Newtonian solutions based on impact theory for hypersonic vehicles 
were calculated using the code presented in the appendix. This code was wTitten 
specifically to be compatible with LAURA. 

Cases 

Four freestream conditions are considered, Mach 4.6, 7, 8.5, and 10, at 5 de- 
grees angle of attack. Complete freestream descriptions are contained in Table 1. 
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Figure 1: Side view of singularity-free surface mesh for Stardust reentry capsule 
forebody. 
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Table 1: Freest ream conditions. 


The Mach 4.6 solution is for a perfect gas, while equilibrium air is used for the 
Mach 7-10 results. 

Three axisymmetric solutions were also obtained, at Mach 7 and 10, equi- 
librium air, and Mach 7, perfect gas. 

All Newtonian results are for perfect gas. 

Results 

Axial-force coefficients are presented in Table 2 for all cases. The reference 
area is 0.5189 m 2 . For Ca , unmodified-Newtonian results with a base-pressure 
correction, as described in the appendix, are tabulated for comparison with the 
LAURA solutions. Excellent, agreement is seen for axial force between the two 
methods for Mach 7 10. At Mach 4.6, 5 degrees angle of attack, the Newtonian 
axial force is 8.5 percent higher than the LAURA result, indicating a breakdown 
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method 

a 

4.6 

7 

Moo 

7 (pg) 

8.5 

10 

LAURA 

0 

— 

1.506 

1.454 

— 

1.515 

Newtonian 

0 

— 

1.51 

— 

— 

1.50 

LAURA 

5 

1.41 

1.477 

— 

1.496 

1.498 

Newtonian 

5 

1.53 

1.50 

— 

1.50 

1.49 


Table 2: Axial-force coefficients. 


of the assumptions behind impact theory at those conditions. 

Looking at the effect of gas model on the Mach-7, axisymmetric solution, 
the perfect-gas Ca is 3.4 percent lower than the equilibrium-air result. 

Normal-force coefficients calculated at 5 degrees angle of attack are plotted 
in Figure 2 versus freestream Mach number. Included with the LAURA results 
are modified-Newtonian solutions. The modified-Newtonian results are nearly 
independent of Mach number at (7/v — 0.04. LAURA predicts the Newtonian 
result at Mach 10, but shows nearly a 50 percent decrease in normal force 
coefficient by Mach 7, with a similar number at Mach 4.6 of (7 jv = 0.021. 

A similar shift in values occurs for the pitching-moment coefficient, Fig- 
ure 3. At Mach 10 the LAURA solution agrees with modified-Newtonian, giving 
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Figure 2: Normal-force coefficient versus Mach number. 
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Figure 3: Pitching-moment coefficient versus Mach number . 


Cm — —0.014. Between Mach 10 and 7 the static stability of the capsule is more 
than halved to -0.0056. Again, the solution at Mach 4.6 has aerodynamics simi- 
lar to Mach 7. This drop in pitching-moment coefficient between Mach 10 and 7 
makes the reentry capsule 60 percent less stable in just 8 seconds of flight time. 
A reduction in static stability for this capsule would likely lead to increased am- 
plitude of angle-of- attack oscillations. The reference length is 0.8128 m and the 
center of gravity is on the centerline 0.285 m back from the nosetip. Pitching 
moments are calculated about the center of gravity. 

Reduction of Static Stability 

Seeking the physical fluid phenomenon causing this abrupt decrease in stabi- 
lizing pitching-moment coefficient, symmetry-plane surface pressure coefficients 
are plotted in Figure 4 for the LAURA Mach 7 and 10 cases. The modified- 
Newtonian solution at Mach 10 is included as well, but the Mach 7 Newtonian 
solution is omitted because it agrees with the Mach 10 result to within 0.5 per- 
cent. 

The trends in the LAURA Mach 10 solution are similar to the results of 
impact theory, with nearly flat profiles along the conical portion of the forebody. 
At Mach 7 the leeside shows more of a recompression than the Mach 10 results, 
though the profile is relatively flat. Leeside pressure coefficients are 5 percent 
higher at Mach 7 than Mach 10. The windside character at the two Mach 
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Figure 4: Symmetry-plane surface pressure coefficients at 5 degrees angle of 
attack. 


numbers differs significantly, with the pressure dropping off significantly toward 
the shoulder at the lower Mach number. This change in character of the surface 
pressure distribution between Mach 10 and 7, an increase in the leeside pressure 
and a dropoff of the windside pressure toward the shoulder, is responsible for 
the reduction in pitching-moment coefficient at these hypersonic speeds. 

Looking at surface pressure distributions, the Mach 7 and 10 LAURA surface- 
pressure coefficients are presented side-by-side in Figure 5. On the conical pro- 
tion of the forebody the Mach 10 contours form predominantly radial lines, 
consistent with impact theory. At Mach 7, however, the surface C P contours 
are highly curved, more consistent with a locally subsonic distribution. The 
patterns seen in Figure 5 suggest a change in the nature of the sonic bubble 
between Mach 7 and 10. 

Figure 6 displays such a change in the nature of the sonic bubble. Plotted 
over a head-on view of the heatshield surface is the location of the sonic line, out- 
side the boundary layer, for the Mach 7 and 10 LAURA solutions. At Mach 10 
the sonic bubble contains the spherical nosepoint and the windward third of the 
forebody cone. By Mach 7, though, the sonic bubble encompasses two-thirds of 
the forebody, allowing significant three-dimensional relieving effects for the fore- 
body flow, resulting in the surface pressure distributions previously compared 
in Figure 5. 
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Summary of Results 

Aerodynamic coefficients are calculated for the Stardust reentry-capsule fore- 
body at Mach 4.6, 7, 8.5, and 10. Axial-force coefficients agree well with 
Newtonian-flow assumptions for Mach 7-10. Normal-force and pitching-moment 
coefficients agree with Newtonian at Mach 10, but reduce sharply in magnitude 
with decreasing Mach number, dropping more than half their values by Mach 7. 

The reduction in static stability by 60 percent is investigated and found to 
correspond with a change in the inviscid sonic-bubble location. The sonic bubble 
grows with increasing Mach number to encompass two-thirds of the forebody 
by Mach 7, at 5 degrees angle of attack. 
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Appendix 

A Modified-Newtonian Surface Pressure 
Calculator Compatible with LAURA [10] 


Abstract 

Modified-Newtonian surface pressure coefficients are calculated and inte- 
grated to obtain lift and drag coefficients. Three-dimensional geometries defined 
as a structured surface mesh are assumed, with freestream pressure imposed on 
shadowed portions of the geometry. Inputs required are the surface mesh, flow 
angle of attack, freestream Mach number, ratio of specific heats, and a nor- 
malizing length and area. Axial/normal and lift/drag force coefficients and the 
pitching moment coefficient are output. A plot file is created of the surface 
pressure coefficient. Compatability with LAURA is emphasized and the source 
code is presented in a User’s Manual format to aid in customization. 


Governing Equations 

Modified-Newtonian surface pressures (see Anderson[ll] §3. 2-3. 5, or Bertin[12] 
§6. 2-6. 3), based on impact theory, are calculated according to, 


Cp 

Cp ,max 


— sin 2 0 = cos 2 cj) 


( 1 ) 


where the pressure coefficient is defined as, 


C P 


P-Poo 


( 2 ) 


Straight-Newtonian is obtained in Eqn. 1 by setting Cp tmax = 2. 6 is the 
inclination angle of the surface to the freestream velocity vector, . <f> is the 
angle between and the unit-normal to the surface, n. <j> and 0 are related 
as 6 = | — 4>. 

Eqn. 1 is evaluated at a point using the relation, 


cos (p = —ey • n (3) 

where ev is the unit vector along V^, so that Voo — T Uo e v. 


Initializations 

The surface is read as a structured mesh in i and j. z, y , and z are the 
nodal coordinates. 
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Figure 7: Orientation of angles for Newtonian flow solver. 


program newtonian 

parameter ( imx = 59, jmx = 30, kmx = 1 ) 

parameter ( itdim = 2 * (imx - 1) * (jmx - 1) ) 

dimension x(imx, jmx, kmx) , y (imx, jmx, kmx) , z( imx, jmx, kmx) 

dimension dxi(imx, jmx) , dyi(imx, jmx) , dzi (imx , jmx) 

dimension dxj (imx, jmx) , dyj (imx, jmx) , dzj (imx, jmx) 

dimension ev(3) , sn(imx , jmx , 3) 

dimension cpcpmax(imx , jmx) , cp(imx,jmx) 

dimension tarea(itdim) , tnormal (itdim, 3) , tcpcpm(itdim) 

dimension tcp(itdim), txyz(itdim,3) 

real mag 

integer diagnos 

diagnos = 0 

From the angle-of-attack, a, obtain the unit vector, ey*, defining the freestream 
velocity direction. In Laura coordinates. 

write (6,*) ’Enter angle-of-attack in degrees’ 

read (5,*) alpha 

pi * acos( -1. ) 

alpha = alpha * pi / 180. 

ev(3) = -cos( alpha ) 

ev(l) = sin( alpha ) 

ev(2) = 0. 

Read the surface grid in single-block Plot3d ASCII 3-D whole format. 

open ( 10, file = ’newton. g’, form = ’formatted’ ) 
read (10,*) idim, jdim, kdim 

if ( idim .gt. imx .or. jdim .gt. jmx .or. kdim .gt. kmx ) then 
write (6,*) ’Increase dimensions: idim, jdim, kdim 


10 


j 


i 



stop 
end if 

read (10,*) ((( x(i,j,k), i=l,idim), j=l,jdim), k=l,kdim), 
$ ((( y Ci , j »k) , i=l , idim) , j=l,jdim), k=l,kdim), 
$ ((( z(i,j,k), i=l,idim), j=l,jdim), k=l,kdim) 


Surface Normals at a Point 


n is computed at each surface grid node by taking the cross product of the 
surface tangents in i and j. 

tl * tj (4) 


n = 


I X tj II 

Surface tangents are computed with second-order discrete derivative approx- 
imations, with no correction for grid stretching. 


At interior nodes, 


At iffiin Of imaxi 


t i = (A Xi, A y u A zi) 

1 


A Xi 


Ari = 


-3a-] + 4ar 2 - £3 


Do (-direction derivatives from min to max. 


do j = 1, jdim 

dxi (1 , j) = ( -3. * x(l,j,l) + 4. * x(2, j ,1) - x(3,j,l) )/2 

dyi(l,j) = ( -3. * y(l, j , 1) + 4. * y(2,j,l) - y(3,j,l) )/2 

dzi(l.j) = ( -3. * z(l,j,l) + 4. * z(2,j,l) - z(3,j,l) )/2 

end do 

do i = 2, idim - 1 
do j = 1, jdim 

dxi(i,j) = ( x(i+l , j , 1) - x(i-l, j ,1) )/2. 

dyi (i > j ) = ( y(i+l,j ,1) - y(i-l,j»D )/2. 

dzi(i,j) = ( z(i+l, j ,1) - z(i-l , j , 1) )/2. 

end do 
end do 

do j = 1, jdim 

dxi(idim,j) = ( 3. * x(idim,j,l) - 4. * x(idim-l , j , 1) 

$ + x(idim-2, j , 1) )/2. 

dyi (idim, j ) = ( 3. * y (idim, j , 1) - 4. * y (idim-1 , j , 1) 

$ + y (idim-2 , j ,1) )/2. 

dzi (idim, j ) = ( 3. * z(idim,j,l) - 4. * z(idim-l , j , 1) 

$ + z(idim-2, j , 1) )/2. 

end do 

Do j -direction derivatives. If f = 1 is a singularity point, set imin ~ 2 and 
assume ty| i= i = tj| i=2 . 
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imin = 1 

do i = imin, idim 

dxj(i,l) = ( -3. * x(i,l,l) + 4. * x(i,2,l) - x(i,3,l) )/2. 

dyj(i,l) = ( -3. * y(i,l,l) + 4. * y(i,2,l) - y(i,3,l) )/2. 

dzj (i, 1) = ( -3. * z(i , 1 , 1) + 4. * z(i ,2 , 1) - z(i,3,l) )/2. 

end do 

do j =2, jdim - 1 
do i = imin, idim 

dxj(i,j) = ( x(i, j+1,1) - x(i,j-l,l) )/2. 
dyj (i , j) = ( y(i, j+1,1) - y(i,j-l,l) )/2. 

dzj (i, j) = ( z(i, j+1,1) - z(i,j-l,l) )/2 . 

end do 
end do 

do i = imin, idim 

dxj (i , jdim) = ( 3. * x(i,jdim,l) - 4. * x(i , jdim-1 , 1) 

$ + x(i , jdim-2 , 1) )/2. 

dyj (i , jdim) * ( 3. * y(i,jdim,l) - 4. * y (i , jdim-1 , 1) 

$ + y (i , jdim-2, 1) )/2. 

dzj(i,jdim) = ( 3. * z(i,jdim,l) - 4. * z(i , jdim-1 , 1) 

$ + z(i, jdim-2, 1) )/2. 

end do 

if ( imin .eq. 2 ) then 

write (6,*) ’handling i=l as pole’ 
do j = 1, jdim 

dxj (1, j) = dxj (2, j) 

dyj(l.j) = dyj ( 2 , j) 

dzj(l.j) = dzj (2, j) 
end do 
end if 

Form the surface normals n. 

do i = 1 , idim 
do j = 1, jdim 

sn(i , j , 1) = dyi(i,j) * dzj(i.j) - dzi(i.j) * dyj(i.j) 

sn(i, j ,2) = dzi(i, j) * dxj(i,j) - dxi(i.j) * dzj(i.j) 

sn(i, j ,3) = dxi(i, j) * dyj(i.j) - dyi(i,j) * dxj(i,j) 

mag = sqrt (sn(i , j , 1) **2 + sn(i,j,2)**2 + sn(i , j ,3) **2) 
sn(i,j,l) = sn(i,j,l) / mag 
sn(i , j ,2) = sn(i , j ,2) / mag 
sn(i,j,3) = sn(i,j,3) / mag 
end do 
end do 
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Cp/Cp. max at a Point 

If the point is shadowed, <f> > j, then Cp = 0. Otherwise apply Eqns. 1 
and 3. 


do i = 1, idim 
do j = 1, jdim 

cosphi = - ( ev(l) * sn(i,j,l) + ev(2) * sn(i,j,2) + 
$ ev(3) * sn(i , j ,3) ) 

cosphi = raax( 0., cosphi ) 
cpcpmax(i,j) = cosphi**2 
end do 
end do 


C 'pjnax 

Compute Cp t max assuming the surface pressure equals the total pressure behind 
a normal shock in Eqn. 2. Using perfect gas assumptions, 



write (6,*) ’Enter Mach number and gamma’ 

read (5,*) amach, gamma 

amach = amach* *2 

gl = gamma - 1 

g2 = gamma + 1 

g3 » gamma 

p02pi = g2**(g2/gl) * amach * 0.5 / 

$ (2. *(2. *g3-gl/ amach ) )**Cl./gl) 

cpmax =2. / g3 / amach * ( p02pi - 1. ) 
if ( cpmax .It. 1. .or. cpmax .gt. 2. ) then 
write (6,*) ’*** Problem with Cpmax = cpmax 
else 

write (6,*) ’Cpmax = cpmax 
end if 

do i = 1 , idim 
do j = 1, jdim 

cp(i,j) = cpcpmax(i,j) * cpmax 
end do 
end do 

Write Tecplot[ 13] data file. 
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open ( 11, file = ’newton.dat’, form = ’formatted’ ) 
write (11,*) ’TITLE = "Modified Newtonian Surface Pressure’, 

$ ’ ratios"’ 

write (11,*) ’VARIABLES = "X", "Y", "Z", "Cp/Cpmax", "Cp"’ 
write (11,*) ’ZONE T - "Surface", I=’,idim, ’ J=’,jdim, 

$ ’ F=P0INT’ 

do j = 1, jdim 

do i = 1 , idim 

write (11,*) x(i , j , 1) , y(i,j,l), z(i,j,l), cpcpmax(i , j) , 

$ cp(i,j) 

end do 
end do 

Integration of Surface Pressures 

Integrate the pressure distribution over the surface to obtain forces and 
moments. To do so, Cp is needed to be defined over surface area elements, not 
just at surface points. 

Surface Areas and Normals 

Triangulate each quadrilateral surface cell (i + l,j), (f+l,j + l), (i,j + l)] 

using a diagonal cut to create triangular cells [(i, j), (i + 1, j), (i, j + 1)] and 
[(i + 1, j 4- 1), (i, j + 1), ( i + 1, j)]. Surface normals are computed using Eqn. 4. 
Surface tangents are computed using first-order differences. 

icell = 0 

do i = imin, idim - 1 
do j = 1, jdim - 1 
icell = icell + 1 

txyz(icell , 1) = ( x(i+l,j,l) + x(i,j,l) + x(i,j+l,l) ) / 3. 

txyz(icell,2) = ( y(i+l,j,l) + y(i,j,l) + y(i,j+l,l) ) / 3. 

txyz(icell , 3) = ( z(i+l,j,l) + z(i,j,l) + z(i,j+l,l) ) / 3. 

tdxi = x(i+l,j,l) - x(i,j,l) 

tdyi = y (i+1 , j , 1) - y(i,j,l) 

tdzi = z(i+l,j,l) - z(i , j , 1) 

tdxj = x(i , j+1 , 1) - x(i , j , 1) 

tdyj = y (i , j + 1 , 1) - y(i,j,l) 

tdzj = z(i , j+1 , 1) - z(i , j , 1) 

tnormal (icell , 1) = tdyi * tdzj - tdzi * tdyj 

tnormal (icell , 2) = tdzi * tdxj - tdxi * tdzj 

tnormal (icell, 3) = tdxi * tdyj - tdyi * tdxj 

tarea(icell) = 0.5 * sqrt( tnormal (icell , 1) **2 + 

$ tnormal (icell ,2) **2 + tnormal (icell ,3) **2 ) 

if ( tarea(icell) .It. 0. ) then 

write (6,*) 'problem: tarea( icell, ' ) = tarea(icell) 
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stop 
end if 

tnorinal(icell, 1) = tnormal (icell, 1) / 2. / tarea(icell) 

tnormal (icell, 2) = tnormal (icell ,2) / 2. / tarea(icell) 

tnormal (icell, 3) = tnormal (icell ,3) / 2. / tarea(icell) 

end do 
end do 

do i = idim - 1, 1, -1 
do j = jdim - 1, 1, -1 
icell = icell + 1 

txyz (icell , 1) = ( x(i,j+l,l) + x(i+l, j+1,1) + x(i+l , j , l))/3 

txyz (icell,2) = ( y(i,j+l,l) + y(i+l,j+l,l) + y (i+1 , j , 1) ) /3 

txyz (icell , 3) = ( z(i,j+l,l) + z(i+l,j+l,l) + z(i+l , j , 1) ) /3 

tdxi = x(i , j+1 , 1) - x(i+l, j+1,1) 

tdyi = y(i, j+1, 1) - y(i+l,j+l,l) 

tdzi = z (i , j+1 , 1) - z(i+l,j+l,l) 

tdxj = x(i+l , j , 1) - x(i+l, j+1,1) 

tdyj = y(i+l,j,l) - y(i+l, j+1,1) 

tdzj = z (i+1 , j , 1) - z(i+l, j+1,1) 

tnormal (icell , 1) = tdyi * tdzj - tdzi * tdyj 

tnormal (icell, 2) = tdzi * tdxj - tdxi * tdzj 

tnormal (icell , 3) = tdxi * tdyj - tdyi * tdxj 

tarea(icell) = 0.5 * sqrt ( tnormal (icell , 1) **2 + 

$ tnormal (icell ,2) **2 + tnormal (icell ,3) **2 ) 

if ( tarea(icell) .It. 0. ) then 

write (6,*) ’problem: tarea( icell, ’) = tarea(icell) 
stop 
end if 

tnormal (icell , 1) = tnormal (icell , 1) / 2. / tarea(icell) 
tnormal (icell , 2) = tnormal (icell ,2) / 2. / tarea(icell) 
tnormal (icell , 3) = tnormal (icell ,3) / 2. / tarea(icell) 
end do 
end do 


Cp/Cp yTnax over an Area 

Compute surface pressure coefficient in a manner analogous to the computation 
at a grid point. 

icellmax = 2 * (idim - 1) * (jdim - 1) 
if ( icellmax .ne. icell ) then 

write (6 , *) ’ problem , icellmax= ’ , icellmax , ’ icell= ’ , icell 
stop 
end if 

if ( icellmax .gt. itdim ) then 

write (6,*) ’problem, icellmax= icellmax, ’ itdim=’ , itdim 
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stop 
end if 

do ic = 1, icellmax 

cosphi = - ( ev(l) * tnormal (ic , 1) + ev(2) * tnormal (ic , 2) 

$ + ev(3) * tnormal(ic ,3) ) 

cosphi s max( 0., cosphi ) 
tcpcpm(ic) = cosphi**2 
tcp(ic) = tcpcpm(ic) * cpmax 

if (diagnos .eq. 2) write (6,*) ’Cell=’, ic,’ Cp=’, tcp(ic) 
end do 


Force Coefficients 

The force coefficients are defined for a resultant force F as, 


Cp hPooVIS 

where S is the normalizing area. Assuming P ^ for the base pressure, the force 
coefficient can be found as, 

-CpA E n E • i F 
S 

Elements e 




where A E is the area of the surface element and i F is the unit vector in the 
direction of F. Using LAURA coordinates. 


write (6,*) ’Enter area S to normalize force coefficients’ 
read (5,*) sarea 
ca = 0 . 
cn = 0. 

do ic = 1, icellmax 

ca = ca + tcp(ic) * tarea(ic) * tnormal (ic, 3) ! Axial 

cn = cn - tcp(ic) * tarea(ic) * tnormal(ic , 1) ! Normal 

end do 

ca = ca / sarea 


cn = cn / sarea 

cd = ca * cos ( alpha ) + cn * sin (alpha ) 
cl = cn * cos ( alpha ) - ca * sin( alpha ) 
write (6,*) ’Force Coefficients:’ 

write (6,*) ’C_A = ’ ™ 

write (6,*) ’C__N = 

write (6,*) ’C_D = 

write (6,*) ’C_L = 


ca 

cn 

cd 

cl 


! Drag 
! Lift 
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Base Pressure Correction 


Base pressures are corrected for the axial force coefficient using the method in 
A PAS of Bonner et a/[ 14]. 


ACU = 


A/ 2 

oo 


0.57 

M 4 

i (X) 


cacor = 1. / amach - 0.57 / amach / amach 
ca = ca + cacor 

cd = ca * cos( alpha ) + cn * sin( alpha ) ! Drag 

cl = cn * cos ( alpha ) - ca * sin( alpha ) ! Lift 

write (6,*) ’Force Coefficients with base pressure correction 

write (6,*) ’C_A = ’ , ca 

write (6,*) ’C_N = ’ , cn 

write (6,*) ’C_D = ’ , cd 

write (6,*) ’C_L = ’ , cl 


Pitching Moment 


The pitching-moment coefficient is defined as, 


Cm = 


M 


hPooVISc 


where c is a normalizing length. The moment coefficient is computed about the 
center of gravity (x ao ,x no ) as, 


Cm = ~ C P A E [(n E * i a )(x nB - x no ) - (n E - i n)(x aE - ar ao )] 


In Laura coordinates, the moment formula becomes, 



Y2 CpAe [(he • i z )(x E - x 0 ) - (n E • i x )(z E ~ «o)] 

E 


write (6,*) ’Enter eg location in LAURA coordinates, (x,z)’ 
read (5,*) xO, zO 

write (6,*) ’Enter normalizing length for moment’ 
read (5,*) alength 
cm = 0. 

do ic = 1, icellmax 

cm = cm + tcp(ic) * tarea(ic) * ( 

$ tnormal(ic,3) * ( txyz(ic,l) - xO ) - 

$ tnormal (ic , 1) * ( txyz(ic,3) - zO ) ) 

end do 

cm = cm / sarea / alength 

write (6,*) ’Moment coefficient, nose up positive:’, cm 


stop 

end 
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